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Euler turbulence has been experimentally observed to relax to a metaequilibrium state 
that does not maximize the Boltzmann entropy, but rather seems to minimize enstrophy. We 
show that a recent generalization of thermodynamics and statistics due to Tsallis is capable 
of explaining this phenomenon in a natural way. The maximization of the generalized 
entropy S1/2 for this system leads to precisely the same profiles predicted by the Restricted 
Minimum Enstrophy theory of Huang and Driscoll. This makes possible the construction of 
a comprehensive thermodynamic description of Euler turbulence. 
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The kinetic and field equations for the drift motion of a pure electron plasma column in a strong magnetic 
field are isomorphic to the equations of motion of a two-dimensional Euler fluiotU. The density of the 
plasma corresponds to the vorticity of the fluid, and the electrostatic potential corresponds to the stream 
function. This observation has made it possible to use pure electron plasmas to study Euler turbulence in the 
laboratoryn. Such experiments have followed the relaxation of Euler turbulence through several identifiable 
stagesu: An initially hollow vorticity profile develops a linear diochotron instability which saturates with the 
creation of long-lived vortex patches. These patches move about for hundreds of diochotron periods, shedding 
filiaments, and eventually mixing and inwardly transporting. This process gives rise to an axisymmetric 
metaequilibrium state, whose density decreases monotonically with radius, which then persists for tens of 
thousands of diochotron periods. The eventual decay of this state is due-only to viscous and three-dimensional 
effects that destroy the idealization of the two-dimensional Euler fluidlm 

The shape of the radial vorticity profile of the metaequilibrium state is an interesting and fundamental 
problem. One would expect that it could be described by a variational principle, but the most natural 
principle of this sort - the maximization of the Boltzmann entropy under the constraints of constant mass, 
energy, and angular momentum - has been found to yield profiles that are substantially flatter than those 
observed in experimentsuO. On the other hand, an alternative variational principle, in which the enstrophy 
(the integral of the-square of the vorticity) is minimized has been found to yield results in excellent agreement 
with experimented. To date, however, there has existed no satisfactory theoretical explanation for this 
unusual variational principle. 

In this paper we show that the failure of the Boltzmann entropy to predict the radial density profile of 
the metaequilibrium state can be understood as but one example of a systemic breakdown of Boltzmann.- 
Gibbs statistics for systems with long-range interactions, long-time memoryQ, or fractal space-time structured. 
Moreover, we show that a recent generalization of statistics and thermodynamics due to Tsallisa is capable of 
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explaining this phenomenon much more naturally: The maximization of the Tsallis entropy S q for this system, 
with q = i, leads to precisely the same profiles predicted by the Restricted Minimum Enstrophy (RME) 
theory of Huang and DriscolE. This observation makes it possible to develop a consistent thermodynamic 
description of such systems, and to associate this phenomenon with a wide body of research on generalized 
statistics and thermodynamics. 

The outline of this paper is as follows: In Section |l| we describe the dynamical equations of the nonneutral 
plasma column (or, equivalently, of the two-dimensional Euler fluid), cast them in Hamiltonian format, and 
present the constants of the motion. We also review the experimental results for this system, and d escr ibe 



previous attempts to explain the metaequilibrium density profile by a variational principle. In Section III 



describe Tsallis' generalization of thermodynamics, and in Section [V we review the application of Tsallis' 
formalism to the problem of stellar polytropes, which are-static solutions to the Poisson-Vlasov equations. 
This problem was first considered by Plastino and PlastinoO, who applied Tsallis' methods to a linear eneijgy 
functional. Noting that the energy functional of the Poisson-Vlasov system is, strictly speaking, quadratidI9, 
we redo this analysis. Finally, in Section ^ we return to the problem of the metaequilibrium state of two- 
dimensional Euler turbulence, and we show that Tsallis' generalized thermodynamics may be used to explain 
the observed density profiles. 



II. NONEUTRAL PLASMA PROFILES AND EULER TURBULENCE 



A. Dynamical Equations 

Strongly magnetized pure-electron plasmas in "Penning .traps" with cylindrical geometry and electrostatic 
axial confinement have been studied for some time nowQ. Such plasmas typically have a gyrofrequency 
that is much greater than the bounce frequency, which in turn is much greater than the drift frequency. 
That being the case, we can average over the gyro and bounce time scales, and describe the system by the 
two-dimensional drift motion of guiding centers, perpendicular to the magnetic field. 

Since the magnetic field is uniform, the dominant drift mechanism is the E x B drift, given by 

E x B 

VExB - C ^ , 

where B is the applied magnetic field, E is the self-consistent electric field, and c is the speed of light. Since 
this drift velocity is independent of particle thermal velocity, it is possible to project out the velocity degrees 
of freedom in phase space, and thereby write a Vlasov equation directly for the guiding-center density n(r, t), 

dn(r,t) _ , 

0= y J ' +v ExB -Vn r,t . 
at 

Writing E = — V$, and adopting dimcnsionlcss units with a magnetic field of unit magnitude, this can be 
written 

= ^ + t. • [V$ x Vn(r. t)\ . (II.l) 

at 

where b is a unit vector in the direction of the magnetic field, and where we have used the vector "triple- 
product" identity. The self-consistent electrostatic potential is then given by the Poisson equation, 

V 2 $(r, t) = 47ren(r,i), 

where — e is the electronic charge. (Henceforth, we set e = 1.) 

If we identify n as the vorticity and $ as the stream function, we note that these equations are isomorphic 
to Euler's equations of inviscid fluid dynamics in two dimensions. Likewise, Dirichlet boundary conditions, 
for which the wall is an equipotential, correspond to the condition that the normal velocity of the Euler 
fluid vanishes at the wall. Indeed, simulations of pure-electron plasma columns provide an important ex- 
perimental tool for the study of two dimensional Euler turbulence. Henceforth, we interchangibly refer to 
the physical embodiment of this dynamical system as a pure-electron plasma column, or as two-dimensional 
Euler turbulence. 
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B. Hamiltonian Structure 



In spite of the fact that we have projected from phase space to configuration space, we note that the 
Vlasov equation, Eq. (II. 1), has a symplectic Hamiltonian form in two dimensions. Specifically, taking the 
magnetic field in the z-direction, so that b = z, the Vlasov equation has the form 

dnjr, t) dn(r,t) &f>(r,t) dn(r,t) d$(r,t) 

~^ = — X dy dy te- = -W'.*). fc ('.*)l. 

where the single-guiding-center Hamiltonian is 

h(r,t) = -$(r,i), 

and we have defined the corresponding Poisson bracket 

da(r,t) db(r,t) da(r, t) db(r, t) 
[a(p,*),fc(r,t)] = -^ g- y g- y sr , 

so that x and y are canonically conjugate variables. The configuration space can thus be regarded as a phase 
space, and the dynamics of the plasma are then a symplectomorphism in configuration space. 

The corresponding Hamiltonian field structure is noncanonical, and Lie-Poisson in formlMx That is, the 
equation of motion is 

^M = {„(r,t),ff[n]}, 

where the field Lie-Poisson bracket of two functionals of n is 

SA SB 



{A,B} = J d 2 r n(r,i) 
and the field Hamiltonian functional is given by 



6n(r, t) ' 5n(r, t) 



H[n] = i /" d 2 r n(r,t)h(r,t) = ~ J d 2 r n(r,t)$(r,t) 

= -\J d 2 r n(r, t) J d 2 r' n(r', t)G(r, r'), (H.2) 

where, in turn, G(r, r') denotes the Green's function of the Poisson problem, 

V 2 G(r,r') = 47n5(r - r'). 

The factor of \ in the Hamiltonian prevents double counting of the energy. 
The Lie-Poisson bracket admits the infinite set of Casimir functionalalJ, 

Z^[n] ee J d 2 riP(n{r,t)), 

where ip is any function of its argument. These Casimir functionals commute with any other functional, 
including the Hamiltonian, and hence they are constants of the motion. We may span the set of Casimir 
functionals with analytic ip by the set 

Zj[n] ee - f d 2 rn 1 {r 1 t), 



J 

indexed by the integers j > 1. The Casimir functional Zi is of special importance; it is called the enstrophy, 
since its analog for the Euler fluid is the integral of the square of the vorticity. 

If we further suppose that the Penning trap is cylindrically symmetric, with a grounded outer wall, then 
the Hamiltonian is invariant under rotation and time translation, so that the angular momentum, 

L W =/dWn(r,*), 

and the total energy, H[n], are also good invariants. 



3 



C. Variational Descriptions of the Metaequilibrium State 



In spite of the elegance of this Hamiltonian structure, both laboratory and numerical experiments indicate 
that some of these theoretical invariants are broken, presumably due to collisional effects that are, of course, 
ignored in a Hamiltonian formulation. i_ . 

Stable equilibria, both axisymmetric and nonaxisymmetric, have been observed for this systemuJ. If the 
plasma is initialized with a hollow density profile, however, the spatial gradients will excite diochotron 
(Kelvin- Helmholtz-like) instabilities on a short time scale, which will, in turn, give rise to much longer-lived 
vortex patches. As these patches move about and collide, they shed filiaments of particles which erode 
the vortex patches further, until a metaequilibrium state with a characteristic profile shape is eventually 
reached. This metaequilibrium state can persist for tens of thousanda.£if diochotron periods, until it is finally 
destroyed by viscous effects which are outside the scope of this paperElu. Here we focus on the metaequilibria 
of initially axisymmetric configurations. 

In the course of the above-described evolution, the total mass Z\ and the angular momentum L are well 
conserved. The energy H is reasonably well conserved. The enstrophy Z 2 tends to decrease in more-or-less 
monotonic fashion, and other Casimir invariants, such as the Boltzmann entropy 



S[n] = - J d 2 r n(r) In [n(r)] , 



are badly broken. For this reason, Z\, L, and H, are often xeferred to as robust or rugged invariants, while 
the Zj with j ' > 2 are termed fragile or dissipated invariantsQ. 

It is tempting to try to derive the shape of the final profile from a variational principle. Most work has 
centered on maximizing the Boltzmann entropy, under the condition that the robust invariants are fixedEE 
Using the Boltzmann entropy, one can demand 

= S(S - aZ l -0H- XL), 

which yields the relationship 

- 1 - In [n(r)] + /3$(r) = a + Xr 2 . (11.3) 

Taking the Laplacian of both sides, we arrive at an equation for the density profile, 

- V 2 [lnn(r)] + 47r/3n(r) = 4A. (II.4) 

Unfortunately, tisM observed metaequilibrium density profiles are significantly more peaked than the solutions 
to this equationtffl. _ . 

Matthaeus and MontgomeryEj have suggested that turbulent relaxation follows a Selective Decay Hy- 
pothesis: The approach to equilibrium is governed by the most slowly decaying fragile invariant. In this 
case, because the enstrophy seemSgto be the most slowly .-decaying of all the fragile invariants, it has been 
proposed that nonneutral plasmasu and Euler turbulenceo tend to minimize enstrophy, rather than max- 
imize the Boltzmann entropy, while still respecting the robust invariants. Indeed, if we replace the above 
variational principle with 

= 5(Z 2 - aZ x - 0H - XL), 
then we are quickly led to the relationship, 

n(r)+/3$(r) =a + Xr 2 . (11.5) 
Taking the Laplacian of both sides, we arrive at the linear Helmholtz equation for the density profile, 

V 2 n(r) + 47r/3n(r) = 4A. (II.6) 
The well behaved cylindrically symmetric solutions to this equation are of the form 
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n{r) = vj {nr) + 



(II.7) 



where Jo is the Bessel function, and the constants /i, v, and k, which have replaced /3, A, and the constant 
of integration, are fully determined by the constrained quantities, Zi, H, and L. 

Because the solutions to the above variational problem often predict a negative density near the wall, 
Huang and Driscoll also introduced a Restricted Minimum Enstrophy (RME) models in which the above 
profile is replaced by the cutoff form 

n(r) = / I [Jo(Kr) - Jo(Kro)] f r ° , (II.8) 

w I for ro < r < r w ' 

where r w is the wall radius (which, in dimensionless units, can be set to unity). The constant ro has replaced 
the constant /i; all three constants are still fully determined by the three constraints. This form is justified 
only by the observation that a negative density near the wall cannot be physical, and that nonmonotonic 
profiles are typically subject to diochotron instabilities. 



Huang an d Dr iscoll then carefully comparcdEl experimental data with the profiles generated by Eqs. (II.4), 



( pi.q ), and ([1.8). They found that the data clearly ruled out the maximum Boltzmann entropy profile of 



Eq. (fiT4V The minimum enstrophy profile of Eq. ([1.6) was much better. Best of all was the RME profile 



of Eq. JH.8[ ). The experimental data was clearly consistent with the Bessel function profiles. 

To date, a completely satisfactory explanation of this tendency to minimize enstrophy, rather than maxi- 
mize entropy, does not exist. In the remainder of this paper, we shall show that this phenomenon is consistent 
with a generalization of thermodynamics and statistical physics recently proposed by TsallisO. Though this 
is still not an explanation per se, it certainly makes possible the association of this phenomenon with a much 
larger - and growing - body of research. 



III. GENERALIZED THERMODYNAMICS 

Tsallisi has proposed a generalization of thermodynamics and statistical physics to describe systems with 
long-range interactions, or with long-time memory. For a system with W microscopic state probabilities 
Pi > 0, that are normalized according to 

w 

! = !>, (ni.9) 

i 

Tsallis bases his formalism upon the following two axioms: 
Axiom 1 The entropy of the system is given by 

1 - V W n q k 

s q = k 1 -^^ = ^- l Y J {p,~pi), 

where k and q are real constants. 



Axiom 2 An experimental measurement of an observable O, whose value in state i is Oi, yields the q- 
expectation value, 

w 

i 
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of the observable O. 

It is to be emphasized that these statements are taken as axioms. As such, their validity is to be decided 
solely by the conclusions to which they lead, and ultimately by comparison with experiment. 
We first note that in the limit as q approaches unity we recover the familiar expressions 

Si = -ky^pjhxpi 

i 

and 

i 

whence we may identify k with Boltzmann's constant, ks- More generally, it has been notedEl that k may 
be g-dependent, and need only coincide with Boltzmann's constant for q = 1; for the purposes of this paper, 
however, we disregard that possibility and henceforth adopt units so that k — ks = 1. In any case, it is 
clear that Tsallis' thermodynamics contain the more orthodox variety as a special case. 

The success of thermodynamics and statistical physics depends crucially upon certain properties of the 
entropy and energy, and much effort has been devoted to showing that many of these axe valid for arbitrary q, 
and to finding appropriate generalizations of the rest. Following Tsallis' presentationE£l, it is straightforward 
to verify the following properties: 

Property 1 The generalized entropy is positive. 

That is, we have S q > 0, where equality holds for pure states (3i : pt = 1) and for q > 0. 
Property 2 The microcanonical ensemble has equiprobability. 



To see this, we extremize the generalized entropy under the constraint of normalized probabilities, Eq. (III. 9) 
Introducing the Lagrange multiplier A, we set 

It follows that 

p -[ *(!-«) 

Pi — 

q 

Since this is independent of i, imposition of the constrant, Eq. ( III.9| ), immediately yields pi = l/W. 
Property 3 The entropy is concave for q > and convex for q < 0. 
This follows immediately from the Hessian matrix, 



d 2 



dpidpj 



f w \ 

Sq - A^p 4 = -^~ 2 %, 



which is clearly negative (positive) definite for q > (q < 0). Thus, the generalized entropy is maximized 
for q > 0, and minimized for q < 0. 

Next, we consider the canonical ensemble. If we define a state energy Si, so that the generalized internal 
energy is given by 



then we can extremize S q under the constraint that probability is conserved and that the energy is fixed. 
We find 



G 



Property 4 The canonical ensemble probability distribution is 

Pi = ^r[l-(l-q)p£i}^ 

where we have defined the generalized partition function 

w 



Z 9 =^[l-(l-g)/3 £i ]T^ 



and where we have defined the inverse temperature, (3 = 1/T. 

We note that, in the limit as q approaches unity, we recover the familiar expressions 

Pi 



Zi 
and 

w 

i 

For 5^1, we note that the absolute value of the energy matters - an additive constant in the energy 
spectrum will produce physical effects. Moreover, we note that, for generic real values of q, the above 
expression for pi breaks down if 1 — (1 — q)f3si < 0. In such a situation, state i is thermally forbidden. For a 
positive energy spectrum that is unbounded above, and assuming that [3 > 0, this will happen for sufficiently 
high Ei if q < 1. Thus, in this situation, the Tsallis distribution has a natural cutoff in energy for q < 1. 

More significantly, we note that 

Property 5 The Leg endre-trans form structure of thermodynamics is invariant for all q. 
To see this, we first note that 



'dp\ l-q 



= u q , 



whence we identify the free energy 



It is then possible to verify that 



and it follows that 




F q — U q — TS Ql 



OSq _ 1 

dU q T' 



As Tsallis points out0, these equations lie at the very heart of thermodynamics, and_the fact that they are 
invariant under q is significant. The grand canonical ensemble has also been treatedEj. 

The most striking and significant differences between the generalized thermodynamics and the more usual 
variety have to do with the extensivity of the state variables. If we partition the microscopic states of the 
system into two disjoint subsets, L = {1, . . . , V} and R = {V + 1, . . . , W}, with respective probabilities 
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PL = ^Pt 

i=l 

and 

w 

i-V+l 

then it is straightforward to verify that 

Property 6 The generalized entropy obeys the following generalization of the Shannon Additivity Property, 
S q (pi,...,pw) =P q L S q (—,...,— ) +p q R S q ( ) +S q (p L ,p R ). 

\PL PL J \ PR PR J 

Alternatively, we can consider the total entropy of two completely independent subsystems, A and B. Since 
the subsystems are independent, the probability that their union AlJ B has subsystem A in state i and 
subsystem B in state j is given by 

p^ B =p?pf. 

After a bit of algebra, we find that 

Property 7 The generalized entropy obeys the following additivity rule 

S q (A UB) = S q (A) + S q (B) + (1 - q)S q (A)S g (B), 

and is thus superadditive (entropy of whole is greater than the sum of its parts) for q < 1 and subadditive for 
q>l. 

Likewise, we find 

Property 8 The generalized expectation value of an observable O obeys the following additivity rule 
O q (A UB) = O q (A) + O q (B) + {l-q) [O q (A)S q (B) + O q (B)S q (A)} . 

Note that, in both cases, extensivity is recovered only when 5=1. 

It is believed - but not proven at the time of this writing - that the Tsallis entropy is the only cpe 
for which all of the above properties hold. Moreover, generalized versions of the Boltzmann H-theoremEJa, 
fluctuation-dissipation theoremo, and Onsager reciprocity theorems exist for all q. The formalism is thus 
an important generalization of most of the principal results of thermodynamics and statistical physics. 

Of course, to verify that this generalization is useful, it is necessary to show that it holds for certain 
physical systems with values of q that are different from unity. In the past two years, much work has been 
done along these lines, and the method has been applied with great benefit-to astrophysical problems such 
as stella r pO |lytropesa, Levy flightscj, the specific heat of the hydrogen atornEa, and numerous other physical 
systemsr 4 rr 5 l. For some of these systems, strict inequalities have been proven, demonstrating that q must be 
different from unity in order to obtain a consistent thermodynamic descriptionEI. 

IV. STELLAR POLYTROPES 

A. Hydrostatic Equilibria 

One of the first problems £a which Tsallis' thermodynamics was applied! was that of stellar polytropes, 
first studied by Lord Kelvirc3, and treated in detail by ChandrasekhaiO. Because stellar polytropes are 
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equilibria of the Poisson-Vlasov equations, they are highly relevant to the present study. Therefore, in 
this section, we shall review the previous application of Tsallis' formalism to this problem by Plastino and 
PlastinoEl. We note that they used a energy functional that was linear in the distribution function, whereas 
the full Poisson-Vlasov energy functional is quadratic. In this treatment, we use the full quadratic functional, 
and compare our analysis to theirs. One of the side benefits of this treatment is that it shows the extension 
of Tsallis' second axiom to observables that are quadratic functionals of the distribution. 
A polytropic process has the equation of state 

P = Kp\ 

where P is the pressure, p is the density, and 7 is a constant that can be related to the specific heats. If 
3>(r) denotes the gravitational potential, then the hydrostatic equilibrium is given by 



= -VP - pV$. 



It follows that 



= V p 



,7-1 



7-1 



We now seek solutions with compact support in domain T>. If we require that the density p vanish on the 
boundary dT>, then we must have the following relationship between p and $, 



P 



7-1 



where is the potential on the boundary. The nonlinear Poisson equation for the gravitational potential 
is then 



V 2 * 



(IV.IO) 



where <I> = $^ ^ — $, and C is a constant. The boundary condition is that f = on dV. This equation 
has, for example, spherically symmetric solutions, corresponding to compact spherical configurations of self- 
gravitating mass, that are called stellar polytropes. 

B. Kinetic Equilibria 

As an alternative to the above hydrodynamic description, we can seek polytropic equilibria of the Vlasov 
equation for the mass distribution function f(z,t) where z = (r, v) coordinatizes the phase space of the 
system. As is well known, the equilibria of the Vlasov equations are functions of the constants of the motion. 
We denote the (negative of the) total energy by 



£{z) ee *(r) - -v\ 

so that a marginally confined particle on dT> with zero velocity has £ = 0, and a confined particle has £ > 0. 
Noting that any function of £ is a solution of the Vlasov equation, we examine solutions of the form 



/ 



Qgn-3/2 for g > o 

for £ < 



where 6 is a constant. The mass density of a spherically symmetric configuration is then given by 



p(r) = / dz' f(z')S(r - r') = / dV /(r, v') 



/23< 



dv Attv 



*(r) 



v 

T 



-3/2 



(IV.ll) 
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The integral gives rise to a beta function which can be expressed in terms of gamma functions to finally yield 

(IV.12) 

Comparing this to Eq. ( ["V.IOD , we can identify 



p(r) = frr)*'*8 I>|W). 



7-1 



7 = 1 + 



C. Variational Principle with Linear Energy Functional 



Note that the stellar polytropes comprise a one-parameter family of equilibria, where the parameter is 
7 (or, equivalently, n). We now examine the question of whether or not these polytropic equilibria are 
thcrmodynamically stable, in the sense that they can be obtained from an entropic variational principle, 
and, if so, for what values of the parameter 7 (or n) this is possible. 

Plastino and Plastino have addressed this questiorJj by extremizing Tsallis' entropy for this problem 

5 «[/] = -^7 j dz [/(*)-/*(*)]- 
under the constraints of fixed mass and energy expectation values, 

M q [f] = [dz f q {z) 



U q [f} = J dz f*( z ){£ + *(r) 



(IV.13) 



In fact, they used Mi and XJ\ in their work, because this was before Tsallis had advanced his second axiom 
about expectation values. This issue was subsequently rectified in a paper by TsallisEa, and we present only 
the corrected version here. 

Introducing Lagrange multipliers, the variational problem 



= 6 (Sg + aMg + PUg) 



yields the equilibrium distribution 



m = q 



1 - (q - l)a - (q - l)/3 ^— + *(r) 
{q [l-(q- 1) (a - <J>(°>) + (q - l)/3£(z)] }^ 



(IV. 14) 



For this to be a power law in £, we select a so that 1 — (q — l)(a — <i>( )) = 0, so we can write f(z) = 
DSW-^iz), where D is a constant. The density measured at a point r is then given by the ^-expectation 
value of the spatial delta function, 

Pq {r) = J dz' f q (z')5(r-r') = D q J dV £^(r,V). 
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We see that this corresponds to Eq. (IV. 11 



) if we identify n — § = q/(l — q), or 



U= 2 



1-9 



As pointed out by Plastino and P lastino El, it is known that n must exceed 1/2 in order to avoid the singularity 
in the gamma function in Eq. (IV. 12), but that values in excess of 5 give rise to unnormalizable mass 
distributions, and are therefore unphysical. This means that q € (-co, |). Thus, stellar polytropes cannot 
be described thermodynamically unless q values less than | are used. 

Finally, note that the aforementioned cutoff of the distribution with energy - a generic feature of Tsallis 
distributions with q < 1 - naturally gives rise to the spatial cutoff of the mass distribution, and hence the 
compact nature of the stellar polytrope. 



D. Variational Principle with Quadratic Energy Functional 



Note that the energy fun ctional in Eq. ( IV. 13 ) was regarded as linear in f q (z) in the above analysis. 
Specifically, in deriving Eq. (IV. 14), we wrote the functional derivative of U q with respect to f(z) as 



SU, 



q - = qf q -\z) 



y + *(*■) 



Strictly speaking, this is not correct because the potential <3>(r) depends on f(z), and we did not account for 
this in the above variation—This is precisely the problem of self-consistency of the field - a crucial feature of 
the Poisson-Vlasov systemEJ. To c orrec t this problem, it is best to write the energy as a quadratic functional 
of the distribution, just as in Eq. (II. 2) - only now g-expectation values should be used throughout. Thus, 



U?[f) 



<l: /''(*)y + i / /"(.:) / </:' /" I ' KV(r. r' ). 



where the superscript Q denotes quadratic, the factor of \ in front of the potential prevents double-counting 
of the energy, and G(r,r') is the Green's function for Poisson's equation which satisfies 

V 2 G(r,r') = 47n5(r-r'). 

Note that the functional derivative of U® with respect to f(z) is now 



where we have defined 



which in turn satisfies 



su9 



« - nfl- 1 



y + * 9 (r) 



$ g (r)^ J dz' /V)G(r,r'), 



V 2 $ g (r) -4^(r) 



(IV.15) 



by 



The variational principle thus results in an equation very similar to Eq. ( IV.14| ), except with £ replaced 



* 9 (r) 



where in turn 
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*,(r) = <&<°> - *,(r). 

The resulting expression for p is then a power law in ty q , rather than in 4". We still conclude that 

3 , q 



n 



2 1-q' 



and so the upper bound on q of | still holds. Note that the nonlinear Poisson equation for the gravitational 
potential, 

is now satisfied by ^> q , rather than by 

This exercise might be dismissed as demonstrating little more than the fact that the potential used in 



Subsection 1VC should be interpreted as the ^-expectation value of the Green's function, rather than as 
the usual one. This objection notwithstanding, the derivation using the quadratic energy functional has the 
following virtues: 

• It more clearly shows that the conclusions reached by this method are valid for the system of particles 
in their own self- consistent field. 

• It demonstrates that Tsallis' second axiom extends in a natural way to quadratic functionals of distri- 
butions. 

• It yields the natural generalization of the potential and the density, and shows that the form of Poisson's 



equation, Eq. (IV.15|), relating them is q- invariant. 



• It is generally more consistent with the flavor and spirit of Tsallis' formalism than previous derivations. 

V. GENERALIZED THERMODYNAMIC DESCRIPTION OF EULER TURBULENCE 

We now return to the problem o f der iving the metaequilibrium profiles of relaxed Euler turbulence. We 



redo the calculation of Subsection II C, using the Tsallis prescriptions for the entropy and the robustly 



conserved quantities. The entropy is thus 

S,[n] = -i I jd\ [n(p)-n«(r)], 
and the constraints, expressed in terms of ^-expectation values, are 

Z q [n] = ^Jd 2 r n q (r) 

H q [n] =-\j d2r « 9 ( r ) / <*V n 9 (r')G(r, r') 

and 

L q [n] = J d 2 r r 2 n q (r). 

Setting 

5 (S q - aZ q - (3H q - XL q ) = 0, 

we find 
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n q (r) — q 
9(9-1) 



•/3$ g (r) = a + Ar 2 , 



(V.16) 



where 



satisfies 



$ ? (r) = /" dV n 9 (r')G(r,r') 



V 2 $ g (r) = 47m 9 (r). 



Applying the Laplacian to Eq. ( V.16| ), we obtain the nonlinear Helmholtz equation, 



_ ^ y2 [™ 1_, ( r )] + 47r/3n«(r) = 4A. 



(V.17) 



(V.18) 



Now, the observed particle density p q (r) - i.e., that which is measured in any experiment - is then the 
^-expectation value of the delta function, 

p,(r) = / d 2 r' n q {r')5(r - r') = n q (r). 



In terms of this, Eqs. (V.16), (V.17), and ( V.1S|) can be written 



Pg " (r) ~ g 
q(q - 1) 



+ /35> 9 (r) = a + Xr 2 



V 2 $ g (r) = 47rp 9 (r) 



1 



9(9-1 



+ 47r/?p g (r) = 4A. 



As q — » 1 , it is seen that this reproduces the results of the maximum Boltzmann entropy relationship, 
Eqs. (II. 5) and (IL4). When q = |, on the other hand, we see that, within trivial redefinitions of the 
Lagrange multipliers and the use o f $1/2 instead of 4>, this reproduces the results obtained by minimizing 



the enstrophy, Eqs. (II. 5) and (II. 6), but for a completely different reason. Moreover, just as in the example 
of the stellar polytrope, the cutoff in density at a finite radius ro appears as a completely natural and generic 
feature of the Tsallis distribution, since q < 1, and does not need ad hoc justification. Thus, we conclude 
that all prior observations that have indicated that the relaxation of two-dimensional Euler turbulence tends 
to follow the RME principle of Huang and Driscoll can now be reinterpreted as rather indicating that it 
maximizes the Tsallis entropy S q for q = \- 

We note in passing that there is an easier way to obtain the result that q = \ . Without going through this 
analysis, we note that the enstrophy itself looks rather like (a linear function of) the Tsallis entropy with 
q = 2. Of course, this is misleading, because it is necessary to use ^-expectation values in the extremization 
process. Nevertheless, TsallisH has shown that one result of not using ^-expectation values in the extrem- 
ization process is to effect the transformation q — > l/q. Hence, we are again led immediately to the result 
q = 5 for this system. 



VI. CONCLUSIONS 



The tendency of a two-dimensional Euler fluid to minimize enstrophy - rather than maximize the Boltz- 
mann entropy - during turbulent relaxation to a metaequilibrium state has resisted theoretical explanation 
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to date. In this work, we have shown that density profiles resulting from the Restricted Minimum Enstrophy 
(RME) theory of Huang and Driscoll also maximize the Tsallis entropy. We have thereby provided an alter- 
native way to understand this phenomenon one in which the density cutoff at finite radius emerges quite 
naturally - and to build a consistent thermodynamical and statistical physical explanation for it. In the 
course of doing this, we have verified Plastino and Plastino's upper bound of | on q for the stellar polytrope 
problem using the full quadratic energy functional for the Poisson-Vlasov system; we have also demonstrated 
the g-invariance of the Poisson equation for these systems. 

While still short of a first-principles explanation of the RME model - it would be nice, for example, to have 
an a priori way of knowing why q should be equal to 1 for this system - the observation that RME is consistent 
with Tsallis statistics does effectively associate it with a large and rapidly expanding body of theory. In 
recent years, generalized thermodynamics has been used to describe numerous, widely disparate physical 
systems, with long-range interactions, long-time (non-Markovian) memory, or fractal space-time structure, 
that have resisted previous attempts at a thermodynamic description. It is hoped that this observation will 
stimulate further research in the use of generalized thermodynamics to describe fluid turbulence. 
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